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0.1 Introduction 

We examine the effect of dynamical tides raised by a companion on a solar-type star. 
In these binaries, gravity or g mode oscillations are excited by the companion in the 
radiative region beneath the convective envelope of the star. They become evanescent in 
the convection zone. 

This is of particular interest in connection with the newly discovered planets, some 
of which are found to orbit around solar-type stars with a period comparable to that of 
the high order g modes of the star. One such example is 51 Pegasi (|^; 0). 

Here, we determine the magnitude of the perturbed velocity induced by the tides at 
the stellar surface. We show that, in the case of 51 Pegasi, this velocity is too small to 
be observed. This result is insensitive to the magnitude of the stellar turbulent viscosity 
assumed and is not affected by the possibility of resonance, which occurs when the 
frequency of the tidal disturbance is close to that of some normal mode of the star. We 
also discuss the orbital evolution and synchronization timescales associated with the tidal 
interaction, a detailed calculation of which will be presented elsewhere. 

0.2 Tidal response to a companion in circular orbit 

In the case of 51 Pegasi, observations indicate that the rotational period of the star is 
between 30 and 40 days (§). Since this is much larger than the orbital period (4.23 
days), we neglect the rotational angular velocity of the primary compared to the orbital 
frequency. Quadrupolar tidal forcing thus occurs through potential perturbations with 
periods which are half the orbital period. 

The tides raised on the star are damped by turbulent friction ( |]l0| ) in the convective 
envelope, and by non-adiabaticity arising from heat transport in the radiative interior 
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(PHl). For the periods of interest, these dissipative mechanisms are predominant away 
from and close to resonance respectively. We then calculate the tidal response in reso- 
nance taking radiative damping alone into account. Away from resonance, we should in 
principle include turbulent dissipation to calculate the tide. However, dissipation in the 
convection zone of solar-type stars is weak enough so that the tidal response away from 
resonance is well approximated by assuming adiabaticity. 

The linearized momentum and mass equations governing the response of the non- 
rotating star to the perturbing potential may be written (|^) 

1 o' 

£^ = _AvP' + ^VP-VvI/T, (1) 

P' = -V • (pO , (2) 

where P is the pressure, p is the density, ^ is the Lagrangian displacement vector and 
the primed quantities are Eulerian perturbations. 

In the adiabatic approximation, the energy equation is AS" = 0, where S is the entropy 
per unit mass and A denotes the Lagrangian perturbation. When non-adiabaticity in 
the radiative core is taken into account, this equation takes on the form: pTd (A5) /dt = 
— V • F', where T is the temperature and F' is the perturbed radiative flux. We suppose 
that close to resonance, the response behaves exactly like a free g mode with very large 
radial wavenumber so that WKB theory can be used together with the local dispersion 
relation to estimate V • F' in this last equation. In addition, AS is related to AP and 
Ap through a standard thermodynamic relation. 

Here, only the companion's dominant tidal term is considered in the perturbing po- 
tential (Q). For a system with a circular orbit, this is given in spherical polar coordinates 
(r, ^, if) by the real part of 

*T (r, e, ifi, t) = fr^ F,!"! (cos 0) e™^'^-"*) (3) 

with n = m = 2, Pji™' being the associated Legendre polynomial. Here uj is the orbital 
angular velocity, D is the orbital separation, and / = —GMp/4D^ , with Mp being the 
mass of the companion. The Lagrangian displacement can then be sought in the form: 

p1'"I(cos6') e*™^'^-'^*). (4) 

We can eliminate P' and p' using the non-radial momentum and energy equations. 
The radial momentum and mass equations then reduce to a pair of ordinary differential 
equations for and S^h. 
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where A = dlnp / dr ~ {dlnP / dr) /T i, Fi being the adiabatic exponent. We have defined A 
such that A = A in the adiabatic approximation, and, when radiative damping is taken 
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into account, A — A/{1 + ie) with e — IGacT'^N^ / {bm^ui^ KpPr^) in the radiative core 
and e = elsewhere. Here a is the Stefan-Boltzmann radiation constant, c is the velocity 
of light, K is the opacity and iV^ = —Ag is the square of the Brunt- Vaisala frequency, g 
being the acceleration due to gravity. 

The solution of this system requires two boundary conditions. At the surface of the star 
we take a free boundary. At r = 0, where equations (|^) and (^) have a regular singularity, 
the boundary condition is that the solutions be regular. 

0.3 Numerical results: tidal response and velocity at 
the surface of the star 

The calculations presented in the previous section are applied to a standard solar model 
(||]). We solve numerically the differential equations ^ and (||) using a shooting method 
to an intermediate fitting point. We define x = r/Rc, where Rc is the outer radius of the 
convective envelope. With this notation, the equations are integrated from Xin = 10~^ 
to Xout = 1.00071256. The radiative core extends from x = to a; ~ 0.7. 

In Figure |l|, we plot the spatial distribution of the real parts of mu^r and mcu^fi for an 
orbital period of Pq — 4.23 d away from resonance and in the adiabatic approximation. 
Away from resonance, the magnitude of the imaginary parts of these quantities is much 
smaller than that of their real parts (which is why the adiabatic approximation can be 
used). Therefore, they represent typical values of the radial and horizontal velocities, 
the maximum values being three and six times larger respectively. Since these quanti- 
ties depend on the perturbing mass through the ratio AIp/{Mp + Mq), they have been 
represented in units of this factor. 

We see from Figure || that a companion orbiting around the star with a period Po = 
4.23 d induces a radial velocity at the stellar surface the maximum of which is between 
10~^ and 6 m/s for Mp between 10~^ and 1 Mq. This is at least one order of magnitude 
smaller than the observed velocity. The period of the tidal oscillation corresponding to 
this orbital period is 2.115 d. For the oscillation to have a period of 4.23 d, the orbital 
period would have to be 8.46 d. The maximum perturbed radial velocity at the surface 
of the star induced by the companion would then be between 2 x 10^"^ and 1 m/s for a 
perturbing mass between 10"'^ and 1 M©. These velocities are at least about 50 times 
smaller than the observed ones. These numbers do not depend on the magnitude of the 
turbulent viscosity assumed. We have also checked that, because of evanescence in the 
convective envelope, they are not affected by the possibility of resonance. 

0.4 Discussion and Conclusion 

The planetary companion interpretation has been questioned recently by the reported 
4.23 d modulation in the line profile of 51 Pegasi (§), and the possibility that this 
modulation may be due to g mode oscillations has been considered ([^). We note that, 
according to our results, such a modulation could not be due to g mode oscillations 
tidally driven by a companion. 

From the calculations presented above, it is also possible to compute the various 
timescales associated with the tidal interaction. For Po = 4.23 d and Mp = 10~^ Mq, it 
is found that if the turbulent viscosity of the star is calibrated such that the calculations 
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can account for the observed circularization rates of main sequence solar-type binaries, 
the tidal orbital evolution, circularization, stellar spin up and convective envelope spin 
up timescales are respectively 143, 24, 131 and 18 Gyr. All of these timescales are 
long compared with the inferred age of 51 Pegasi If the companion is a low- 

mass star of O.lAf©, as has been recently suggested, these numbers drop to 1.7, 0.25, 
0.016 and 0.0022 Gyr respectively. We then expect the primary star to be synchronized 
with the orbit, in which case exchange of angular momentum is no longer taking place. 
Synchronization is actually expected if the mass of the companion is larger than about 
10 Jupiter masses. We note that since the orbital decay timescale is larger than the 
synchronization one, tidal interaction stops before the companion has plunged into the 
central star. 

If a simple estimate of the turbulent viscosity based on the usual mixing length theory 
is used, all these timescales have to multiplied by ~ 50. In that case, synchronization is 
expected if the mass of the companion is larger than about 70 Jupiter masses. 
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